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Abstract 

In a previous paper we have set up the Wheeler-DeWitt equation 
which describes the quantum general relativistic collapse of a spherical 
dust cloud. In the present paper we specialize this equation to the case 
of matter perturbations around a black hole, and show that in the 
WKB approximation, the wave-functional describes an eternal black 
hole in equilibrium with a thermal bath at Hawking temperature. 
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1 Introduction 



Quantum gravitational effects are expected to modify the nature of singu- 
larities that arise as the end state of the classical gravitational collapse of a 
compact object. A concrete analytical model of classical spherical collapse is 
the Lemaitre-Tolman-Bondi (LTB) dust solution of Einstein equations, which 
shows that the singularity forming in the collapse could be either naked or 
covered, depending on the choice of initial conditions [1, 2]. 

Treating this dust collapse model as a classical background, one can quan- 
tize a massless scalar field on this space-time using standard techniques 
[3, 4, 5]. When the classical collapse ends in a black hole, the quantiza- 
tion of the scalar field yields the emission of Hawking radiation from the 
black hole, as expected. However, a strikingly different result is obtained 
when the scalar field is quantized on a classical background which ends in a 
naked singularity. It turns out that during the period of validity of the semi- 
classical approximation (curvatures should be less than Planck scale), the 
collapsing cloud emits only about one Planck unit of energy [6]. Moreover, 
because the back-reaction does not become important so long as gravity can 
be treated classically, it follows that the future evolution of the star is gov- 
erned by quantum gravitational effects and it is impossible to say, from the 
semi-classical approximation, whether the star radiates away its energy on a 
short time scale or settles down into a black hole state. This is completely 
different from the black hole case where essentially the entire star evaporates 
via Hawking radiation during the semi-classical phase. 

A full understanding of the gravitational collapse, both in the naked and 
in the covered case, requires the application of quantum gravity. Given our 
limited understanding of quantum gravity at present, perhaps the theory 
which is currently most suited for addressing dynamical quantum collapse 
and questions regarding its final state is canonical quantum general relativ- 
ity. Although limited in its ultimate scope as a theory of quantum gravity, 
the canonical theory can meaningfully address the issue of singularities in 
minisuperspace and midisuperspace models, so long as one can tackle ques- 
tions relating to operator ordering and regularization, and provide a notion 
of time evolution in the theory. 

The minisuperspace quantization of a collapsing null dust shell was ana- 
lyzed in [7] where it was shown how the classical singularity can be avoided 
in the quantum theory. In this model, the avoidance is a direct consequence 



2 



of the unitary time evolution - since the wave-function vanishes at r = 
for early times, it does so at any time. As a consequence, an ingoing quan- 
tum shell develops into a superposition of ingoing and outgoing shell. Such 
a scenario has interesting physical features that cannot be seen in a semi- 
classical approximation: no event horizon can form, and there is neither an 
information loss nor a naked singularity. 

A midisuperspace program to study the canonical quantum dynamics of 
the LTB dust collapse has been developed during the last two years [8], fol- 
lowing earlier pioneering work by Kastrup and Thiemann [9] and Kuchaf 
[10] on the quantization of the Schwarzschild geometry. Here one sets up a 
canonical description of the collapse, using the dust proper time, the area ra- 
dius and the mass function of the cloud as canonical configuration variables. 
The evolution is recorded by the dust proper time. One then develops the 
quantization via the momentum constraint and the Wheeler-DeWitt equa- 
tion, to which a solution for a general mass function has been found using 
an ad hoc delta-function regularization used by DeWitt [11]. We show below 
that this regularization scheme is equivalent to the WKB approximation. 
(A similar analysis has also been carried out for null dust [12] although the 
constraints that are obtained in this case are linear and there is no need for 
regularization.) 

The Schwarzschild black hole can be viewed as a LTB model with a 
constant mass function. We have applied [13] the program described above 
to this particular case and shown how the horizon area quantization and 
entropy of the eternal black hole can be understood in terms of quantized 
shells of matter that are trapped inside the horizon. 

The purpose of the present paper is to show that the WKB solution is 
able also to describe an eternal black hole in equilibrium with a thermal bath 
at Hawking temperature, if the mass function is chosen appropriately. We 
argue that Hawking radiation may be understood as a combination of the 
WKB and the Born-Oppenheimer approximations on the full quantum wave- 
functional. It reinforces our belief in the overall consistency of the program 
and, in particular, suggests that our proposed choice of operator ordering and 
our definition of the inner product on the Hilbert space of wave-functionals 
capture features of the full theory. As we will see below, the definition of the 
inner product enters crucially in the calculation of the Hawking radiation. 

In Section 2, we briefly recall key results from our previous paper on 
quantum dust collapse [8] , leading up to the WKB solution of the Wheeler- 
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DeWitt equation. In Section 3, this solution is specialized to the case of 
matter around a black hole, and shown to describe Hawking radiation. 



2 Canonical quantization of dust collapse 

The spherical gravitational collapse of a dust cloud having energy density 
e(r, p) in an asymptotically flat space-time is described in comoving coordi- 
nates (r, p,0,4>) by the LTB metric 

ds-" = dr^ - ^^^dp^ - R\p, T)dn\ (1) 

1 + Jyp) 

and the Einstein equations 

<r, P) = ^^ R*(r, p) = ± V/ + F/R . (2) 

Here, F{p) is twice the mass to the interior of the coordinate p and -R(r, p) is 
the area radius of the shell labeled p at the dust proper time r. A tilde and 
an asterisk denote partial derivatives with respect to p and r, respectively. 
(Throughout this paper, the gravitational constant is set equal to one.) 

The canonical dynamics of the collapsing cloud is described by embedding 
the spherically symmetric ADM 4-metric, 

ds^ = - L'idr - N'^dtf - RHn\ (3) 

in the LTB space-time (1), and by casting the action for the Einstein-dust 
system in canonical form. The phase space of non-rotating dust is described 
by the dust proper time, r, and its conjugate momentum, whereas the grav- 
itational phase space consists of the configuration space variables (i?, L) and 
their conjugate momenta. Using a version of the canonical transformation de- 
veloped by Kuchaf [10], the configuration variable L is replaced by a new vari- 
able F (the mass function). In terms of the new chart (r, R, F, P^-, Pji, Pp), 
the momentum and the Hamiltonian constraints read [8] 

Hr = r'Pr + R'Pr + F'Pf « , (4) 
H^P^ + J^P'^-^^0. (5) 
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Here, JF = 1 — F/R. The Hamiltonian constraint shows that on the effective 
configuration space (r, i?), the DeWitt super-metric is just diag(l, This 
is a fiat metric, therefore a redefinition of the area coordinate according to 



R* 



I 



dR 



(6) 



where the positive sign refers to the region exterior to the horizon {R > F) 
and the negative sign to the region interior to the horizon {R < F), brings 
the super-metric to manifestly flat form. In terms of the momentum P^, 
conjugate to R^, the Hamiltonian constraint reads 



H 



p/2 

p-i J^p-i 



. 



(7) 



Quantization is implemented by raising the momenta to operator status and 
requiring the physical state, ^[r, i?*, F], to be annihilated by the constraints. 
In this way, the time evolution of ^[r, R^, F\ is determined by the Hamilto- 
nian constraint, 

p 52 p,2- 

+ 



5r2 ^ 5RI 



AT 



^{r,R,F) = 



(8) 



(where the positive sign before the second term refers to the region outside 
the horizon, and the negative sign to the region inside), while invariance 
under spatial diffeomorphisms is implemented by the momentum constraint. 



^'5t^^'*5R^ 



+ F'± 
SF 



(9) 



In the region exterior to the horizon, Eq. (8) is no longer hyperbolic, in 
contrast to the Wheeler-DeWitt equation on the original configuration space. 
The reason lies in the canonical transformations performed, which lead to a 
new effective configuration space. 

To complete the quantum theory, one must define an inner product on 
the Hilbert space of wave-functionals. In [8], we defined it in a natural way 
by exploiting the fact that the DeWitt super-metric is manifestly flat in the 
configuration space (r, i?*). 



(10) 
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Note that this inner product is defined on a r = constant hypersurface. We 
emphasize that this inner product is in general r-dcpendent. The reason is 
that the Wheeler-DeWitt equation preserves a Klein-Gordon type of inner 
product, not a Schrodinger-type of product [11]. However, as has been shown 
in [14], the Schrodinger inner product is approximately conserved in the 
highest orders of a scmiclassical approximation. Since wc shall deal here 
with WKB states only (Sec. 3), quantum-gravitational correction terms to 
this conservation do not play any role here. Equations (8)-(10) clearly imply 
a specific choice, albeit a natural one, of operator ordering. We will see below 
that this choice is sufficient to reproduce the Hawking effect. 

The momentum constraint is obeyed by any functional that is a spatial 
scalar and, in particular, by the functional 



provided that W has no explicit dependence on the radial label coordinate, 
r. Our choice, while not unique, is dictated by the knowledge that F' is the 
proper energy density of the collapsing cloud. When (11) is substituted in 
the Wheeler-DeWitt equation (8) one finds, on using DeWitt's 5— function 
regularization (5(0) = = (5(")(0) V n e N), that W obeys 



We emphasize again that this regularization prescription is at this stage 
completely ad hoc, and could only be justified from an understanding of the 
full theory. It is even imaginable that, analogous to string theory, Schwinger 
terms may arise in the commutation relations of the constraints that could 
forbid the implementation of the Wheeler-DeWitt equation [15] . Fortunately, 
however, for our present purpose of recovering Hawking radiation it is not 
necessary to resolve this issue. 

Equation (12) yields the solution 




(11) 




(12) 




(13) 
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outside the horizon, and 



/R- /Ftanh-S/- 




(14) 



inside. 



3 Origin of Hawking radiation 

We begin by noting that the quantum constraint (12), which has been writ- 
ten using DeWitt's delta-function regularization, is the same equation as one 
would get by writing the Wheeler-DeWitt equation (8) in the highest-order 
WKB approximation. The reason is that this prescription effectively sup- 
presses the (divergent) WKB prefactor. To show this, let us expand the 
wave-functional ^ of (8) (with Ti being re-inserted) in a power-series in ^, 



^ = exp{iS/h) = exp 



1 r'^ 

— j^ drF'{r)S{T{r),R{r),F{r)) 



(15) 



cS(r, R, F)^So + hSi + h''S2 + ... ■ (16) 

Substituting this expansion in (8) and retaining only the leading order, h- 
independent, terms gives 



dr 



± 



'dSo 



1 



* = 



(17) 



Comparison with (12) shows that the WKB solution is the same as one would 
get by doing the delta-function regularization in the original Wheeler-DeWitt 
equation (after the identification W = iSo). (For a similar discussion of 
WKB states for the Schwarzschild black hole see [16] and for two-dimensional 
dilaton gravity see [17]). 

Wc will now show that the wave- functional (11), along with the solutions 
(13) and (14), yields Hawking radiation when it is applied to a matter distri- 
bution that is appropriate to a massive black hole surrounded by dust whose 
total energy is small compared with the mass of the black hole. For this 
purpose, let us assume that the mass function F{r) is of the form 



F(r) = 2M9{r) + f{r) 



(18) 
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where 6{r) is the Heaviside step-function, and /(r) (not to be confused 
with /(p) occurring in (1)) is differentiable, representing a dust distribu- 
tion with /(r)/2M <^ 1. This mass function is interpreted as the presence of 
a Schwarzschild black hole of mass M at the origin, and /(r) is a dust matter 
perturbation on the black hole, which, as we now show, can be related to 
Hawking radiation. In a sense, it plays the role of the quantum field used in 
standard derivations of Hawking radiation. 

Inserting this mass function in the wave-functional (11) gives (setting 
h — 1 again) 



"1 f'^ 

*[t, R, F] = exp - / dr[2M6(r) + /'(r)] >V(r(r), R(r), F(r)) 
L2 Jo 



(19) 



or 



" 1 POO 

^[t, R, F] = exp [MWoir, R, F)] x exp - / drf(r)W^ Mr), R(r), F(r)) 

L2 Jo 

(20) 

where Wo{t,R,F) = W{t{0), R{0), F{0)). The first exponent on the right- 
hand side is the WKB wave-functional representing the black hole at the 
origin, as shown in [13]. The second term, up to order /(r), represents a 
matter distribution that propagates in this background, if we take F{r) ^ 2M 
in W-^. Thus we have 

^ = ^bhX^/, (21) 



where 



~ exp 



exp[M>Vo(T,i?,F)] , 
1 



I POO 

- drf{r)Wf{r{r),R{r),F{r)) 
I Jo 



The solution is known from (11), (13) and (14) above and given by 



In 



R+V2M' 
VR-V2M 



(22) 



outside the horizon, and 
W/ = -i 



T±2V2M(VR-^ln 



R + V2M' 
\f2M-\fR 



(23) 
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inside the horizon. 

We would hke to rewrite the expressions for W-^ in terms of the Kilhng 
time T. For the Schwarzschild background being considered for the distribu- 
tion /(r) and for contracting clouds, we have the following relation between 
the proper time and the Killing time (see e.g. [18]), 



= T + 



V2M j 



dR- 



R-2M 



T + 2V2M 



R+V2M \ 

Vr-V2m) 



(24) 
(25) 



Thus, in terms of T, there are two possibilities for W/ut, 



out 



T + 4 V2M [y/R - ^ In 



VR-V2M 



(26) 



The wave-functional that corresponds to an inf ailing wave at T — > — cx3 and 
R ^ CO is the one with W/ut given by the second of the above. We will 
therefore concentrate on 



"iff — exp 



; J drf\r) + 4V2M {^\fR - ^^^^ In 



(27) 

Defining Z = AV2MR we find that as — > oo this wave-functional ap- 
proaches 



f ~ exp 



-t J drf'{r)[T + Z] 



(28) 



which undergoes rapid oscillations except when T — > — oo, that is on T~ . 
When T — > oo, in order for the phase to be not large we see that it! — > 2M, 
i.e. Z — > 8M. In this limit we have 



= exp l^-i J drf{r) (t + 



Z-4Mln 

Z - 



exp |-i drf{r) (t + 4Mln 



Z + %M 
Z-8M 



16M 



(29) 



This is similar to what happens in the geometric optics approximation. The 
simple looking phase on X~ has scattered through the geometry to turn into 
the complicated looking phase on X"*" near the horizon. 
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Equation (28) represents infalling waves. We can think of it as a product 
over plane waves, one at each label r, as follows: 

= JJe-^'^MI^W+^Wl . (30) 

r 

This should represent a complete set of infalling modes at each label r, if we 
think of the a;(r) as the frequency of the modes. In other words, we allow all 
possible u;{r) — A/(r). 

A complete set of outgoing modes on X"*" would likewise be given by the 
functional 

^ n e^*'^^'')!^^'')"-^^'')! . (31) 

r 

This is because the transformation from the dust proper time to the Killing 
time is now obtained by matching an expanding (rather than contracting) 
dust cloud to a Schwarzschild exterior, for which one gets, instead of (25), 
the relation 



T -2V2M 



2 \^/R-^/2M 



(33) 



It is then easily shown that as — > 00, the asymptotic form of the wave- 
functional is as in (31). 

Now we ask the question: what is the projection of our solution (29) on 
the negative frequency modes of the outgoing basis on X"*". For this purpose 
we must consider the inner product of states on a hypersurface of constant 
Schwarzschild time T. Thus we must transform from the Euclidean fiat 
metric on the {R^ , r) plane of Eqn. (8) to the metric in the (i?, T) coordinates. 
Using the relations (6) and (25) we get 

(dRA'' ( drV ( R 

9rr = 9R.R* + 7777 9 



^dR J ** \dRj ^ \R-2M 

(Note that in (8) the positive sign holds outside the horizon, so that Qr^r^ — 
+ 1.) The required inner product on a constant T hypersurface then is 

(*-|*;) = n / V9^dR{r)^+{Z{r),T{r),u;(r))^^(Z{r),T(r),f{r)) 
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= 11/ dR^-^n{Z{rmr)Mrm{Z{r),nr),f{r)) 

n / dZ^-^^^{Z,T,co)^^{Z,TJ) . (34) 



r 



This projection represents the negative frequency modes present in the solu- 
tion. We are interested in | /^+|^^\ p because these are the analogs of the 



Bogoliubov coefficients, a;)| . If we think of a;) as 

(3{f,u) = Yl(3{f{r),u{r)), (35) 



we have 



P{f{r),u;{r))^2M [ 



Ism Z -8M 

2M 



-e e 



r dZ {Z - %M)-^-^'^^f e"^^ . (36) 



(16M)-^^^^^ hM 
Substituting u — Z — 8M and integrating we find 

where a = 4MA/. This gives 



A/ 



gSTTMA/ _ I 



(37) 



This is interpreted as the eternal black hole being in equilibrium with a 

thermal bath at the Hawking temperature (SvrM)"^. Our derivation provides 
a functional Schrodinger picture for dust Hawking radiation, consistent with 
the WKB wave-functional which solves the Wheeler-DeWitt equation. 



11 



4 Concluding remarks 



In this paper we have obtained a derivation of Hawking radiation for dust 

matter, starting from the WKB wave-functional which satisfies the Wheeler- 
DeWitt equation for quantum spherical dust collapse. The fact that such a 
derivation could be found should be treated as support for the validity of the 
inner product defined in Equation (10). 

A functional description of Hawking radiation can also be given within 
a Born-Oppenheimer type of approximation to quantum gravity. Instead of 
our wave function one has there a Gaussian quantum state for a quantum 
field. Evolving this state through the background of an object collapsing to a 
black hole, one finds that it encodes information about the Hawking radiation 
similar to our [19]. 

The present analysis also suggests that an exact treatment of the Wheeler- 
DeWitt equation (8) which goes beyond the WKB approximation should 
yield corrections to Hawking radiation, and provide a better understanding 
of the end state of gravitational collapse. This could perhaps be done along 
the general lines presented in [14] in which corrections to the semiclassical 
limit have been calculated from the Wheeler- De Witt equation. These issues 
are at present under investigation. 
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